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ABSTRACT 

We report recent ground-based photometry of the transiting super-Earth exoplanet GJ 1214b at sev- 
eral wavelengths, including the infrared near 1.25 /an (J-band). We observed a J-band transit with the 
FLAMINGOS infrared imager and the 2.1-meter telescope on Kitt, Peak, and we observed several optical 
transits using a 0.5-meter telescope on Kitt Peak and the 0.36-meter Universidad de Monterrey Observa- 
tory telescope. Our high-precision J-band observations exploit the brightness of the M-dwarf host star at 
this infrared wavelength as compared to the optical, as well as being significantly less affected by stellar 
activity and limb darkening. We fit the J-band transit to obtain an independent determination of the 
planetary and stellar radii. Our radius for the planet (2.61!To hi excellent agreement with the 

discovery value reported by Charbonneau et al. based on optical data. We demonstrate that the plane- 
tary radius is insensitive to degeneracies in the fitting process. We use all of our observations to improve 
the transit ephemeris, finding P = 1.5804043 ±0.0000005 days, and To - 2454964.94300 ± 0.00006 BJD. 

Subject headings: stars: planetary systems - transits - techniques: photometric 


1. INTRODUCTION 

A current frontier in extrasolar planet research is the 
detection and atmospheric characterization of super-Earth 
exoplanets orbiting nearby M-dwarf stars (Charbonneau k 
Deming 2007). The recent detection of GJ 1214b (Char- 
bonneau et al. 2009) is a significant step toward this fron- 
tier. GJ 1214b was found by the MEarth survey of nearby 
M-dwarfs (Nutzman & Charbonneau 2008), and it is the 
first super-Earth exoplanet discovered by a ground-based 
photometric method. The radius of GJ 1214b derived from 
the discovery observations (2.68.R©, see Table 1) indicates 
that the planet has a low molecular weight composition 
such as an ocean planet, but alternative compositions 
remain possible (Miller-Rlcci & Fortney 2010; Rogers & 
Seager 2010). Inferences concerning the composition of 
GJ 1214b depend on a securely measured radius. In this 
regard it is valuable to measure the transit in the infrared 
(IR). IR transit observations are less sensitive to stellar 
activity and limb-darkening, compared with optical wave- 
lengths. Reduced sensitivity to these effects helps establish 
a more robust planetary radius. GJ 1214 is also the tar- 
get of several near-term observational programs using both 
HST and Spitzer, so recent ground-based transit observa- 
tions can benefit the space-borne programs by improving 
the precision of the transit ephemeris. 

In this Letter we present results from high-precision pho- 
tometry of a transit of GJ 1214b observed in the J-band 
(1.25 /un), and several less precise transits in the optical. 
Our IR transit data provide an independent check on the 
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planetary radius derived by Charbonneau et al. (2009) 
from the optical discovery observations, and the totality 
of our recent transit, data allow us to significantly improve 
the transit ephemeris. 

2. OBSERVATIONS AND PHOTOMETRY 
2.1. Observations 

We observed transits of GJ1214b using the Universi- 
dad de Monterrey Observatory (UDEM) telescope on UT 
dates April 29, June 6 k 17, 2010. UDEM is a small pri- 
vate college observatory having Minor Planet Center Code 
720, and is located at 689 meters altitude in the suburbs 
of Monterrey, Mexico. The UDEM data were acquired 
using an / c -band filter on the 0.36-meter reflector, with a 
1280x1 024- pixel CCD camera at 1.0 arc-sec per pixel scale. 
We observed an additional transit on UT May 29, 2010 us- 
ing two telescopes on Kitt Peak. We used the 2.1-meter 
reflector with the FLAMINGOS 2048x2048- pixel infrared 
imager (Elston 1998), and a J-band (1.25 /mi) filter, at 0.6 
arc-sec per pixel scale. Following the conclusion of nightly 
public programs, we have access to the 0.5-meter telescope 
at the Kitt Peak Visitor Center (VCT). Simultaneous with 
the 2.1-meter J-band observations, we observed the May 
29 transit using the VCT and a 3072x2048 CCD camera at 
0.45 arcsec per pixel. The VCT observations were planned 
to use the z-band filter, but due to a filter wheel error these 
data were actually acquired in a B filter. Observations at 
all three telescopes used a defocus to improve the pho- 
tometric precision, and all used off-axis guiding to main- 
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tain pointing stability. Exposure times were 30-seconds at 
the 2.1-meter, 60-seconds at the VCT, and 120-seconds at 
UDEM. All of the optical CCD exposures were binned 2x2 
to facilitate rapid readout. 

Flat-field observations were acquired at all three obser- 
vatories using either twilight sky (UDEM and VCT) or 
a series of night-skv FLAMINGOS exposures, incorporat- 
ing pointing offsets to allow removal of stars via a median 
filter. 

2.2. Photometry 

Subsequent to dark current subtraction and division by 
a flat-field frame, we performed aperture photometry on 
the target star and comparison stars. The 20-arcmin field 
of FLAMINGOS provided 8 comparison stars of compa- 
rable IE brightness to G.J 1214. We used 7-pixel circular 
aperture radii to measure the stars. We measured the sky 
background surrounding each star using an annulus with 
an inner radius of 9-pixels and an outer radius of 44-pixels. 
Normalizing GJ 1214 to the the comparison stars yielded 
a transit light curve with an observed scatter that varied 
from 0.002 to 0.0015 as a function of time, due to the de- 
creasing airmass during the observations. We found that 
the best results were obtained by averaging the ratios of 
GJ 1214 to each comparison star. This produced smaller 
scatter than the method of ratioing GJ 1214 to the sum of 
all the comparison stars. We estimated an error for each 
GJ 1214 photometric point as the standard deviation of 
the ratio to the individual comparison stars, divided by 
the square root of their number (error of the mean). The 
scatter in the photometry is approximately 65% greater 
than these estimates, so we increased our error estimates 
by this factor to facilitate more accurate y 2 analyses (see 
below). 

After normalizing to the comparison stars, the GJ 1214 
data clearly reveal the transit but with low-order curvature 
in the out of transit baseline. It is probable that this base- 
line curvature is caused by differences in the effective wave- 
lengths of the J-bandpass as a function of stellar color, in 
combination with the wavelength-variation in telluric wa- 
ter vapor absorption. We account for this baseline using 
a 4th-order polynominal fit, simultaneous with the transit 
fit (see See. 3). Photometry of the UDEM and VCT data 
used similar procedures as for the FLAMINGOS data, ex- 
cept that the baseline was modeled as a linear function 
of airmass. Figure 1 shows the 5 transits that we have 
observed. Figure 2 replots the J-band transit, showing 
both the raw photometry and baseline-removed photome- 
try. and overplots a theoretical transit curve based on the 
(Charbonneau et al. 2009) system parameters. 

3. DERIVATION OF STELLAR AND PLANETARY RADII 

Our J-band transit data have greater photometric preci- 
sion than our optical data (see Figure I). Moreover, stellar 
limb darkening and activity signatures are reduced at IE 
wavelengths compared to the optical. Therefore we use 
only our J-band data to fit for the exoplanet and stellar 
radius, but we use all of our transit data to update the 
transit ephemeris. 

Given an estimate of the stellar mass, three unknowns 
(stellar and planetary radii, and impact parameter) follow 
from three constraints (transit depth and duration, and 


ingress/ egress time). The solution is only weakly affected 
by errors in the adopted stellar mass (Brown et. al. 2001). 
Given also that the IR limb darkening is weak, we fix both 
the stellar mass and limb darkening in our fitting for radii. 
We adopt the stellar mass (0.157 solar masses) given by 
Charbonneau et al. (2009), and the model-atrnosphere- 
based square-root limb darkening coefficients for M-dwarf 
stars derived by Claret (1998). 

We calculate transit curves numerically, using a tile-the- 
star procedure with adaptive-mesh tiles to maximize the 
precision. The transit curve that results from the Char- 
bonneau et al. (2009) system parameters, and the Claret 
(1998) limb- darkening, is illustrated in the lower panel of 
Figure 2. For Figure 2. we overlaid this transit curve on 
our data by first removing the baseline from the data. We 
then varied the central transit time to find the best fit, 
not changing the depth or shape of the curve. It pro- 
vides an excellent fit, note particularly that it accounts for 
the sharp ingress /egress portions of the data. (Figure 2). 
This agreement immediately suggests that the Charbon- 
neau et al. (2009) radii are not signficantly affected by 
limb- darkening errors, or stellar activity. Note also that 
the transit curve still contains significant curvature, un- 
like transit curves at far-IR wavelengths (Richardson et 
al. 2006), where stellar limb darkening is virtually zero. 
The linear limb darkening coefficient for this M-dwarf star 
is about a factor of two lower in J-band than in I-band 
(Claret 1998). 

Given that the Charbonneau et al. (2009) system pa- 
rameters are very close to best-fit values for our J-band 
data, we identify a refined solution as follows. We gener- 
ate a 2-D grid of transit curves centered on the Charbon- 
neau et al. (2009) values, as a function of the stellar radius 
and impact parameter. At each grid point, we find the 
best corresponding exoplanet radius by scaling the depth 
of the transit curve and fitting it to the raw photometry 
simultaneously with the baseline polynominal, using linear 
regression. (Scaling the depth of the transit curve relies 
on the excellent approximation that the transit depth is 
primarily dependent on the ratio of planetary-to-stellar ra- 
dius, and not on the impact parameter.) We also shift each 
transit, curve in time for the best fit. At each grid point 
we calculate the y 2 °f the best fit, and the location of the 
minimum y 2 gives our global best-fit stellar radius, impact 
parameter, and planetary radius. We determine errors on 
these quantities via the range of values contained within a 
X 2 contour of a given significance. The results are given In 
Table 1, except for the best-fit central transit time - that 
is included with the optical transit times in Table 2. 

To fit our four optical transits, we use the Claret (1998) 
limb darkening coefficients to generate synthetic transit 
curves at the optical wavelengths. We shift these in time, 
without varying their depth, to find the best central time 
for each transit. Those times are listed in Table 2. We 
combine our Table 2 transit times with the individual tran- 
sits given by Charbonneau et al. (2009) to derive an up- 
dated ephemeris, listed in Table 1. 

4. RESULTS AND DISCUSSION 

Figure 3 shows our contours of y 2 for the fit to our J~ 
band transit, converted to standard deviations, and shown 
in the plane defined by stellar radius and impact param- 
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efcer. Our x 2 minimum (square on Figure 3) occurs at a 
smaller stellar radius and smaller impact parameter than 
found by Charbonneau et al, (2009). We do not regard 
this difference as physically significant because: 1) it can 
reflect a degeneracy in the fitting process, 2) the difference 
is not much greater than the Charbonneau et al. (2009) 
error range in combination with our errors - the Charbon- 
neau et al. (2009) values are within our 2a contour, and 
3) our radius for the planet is in excellent agreement with 
Charbonneau et al. (2009). 

The contours in Figure 3 trace a parabolic locus wherein 
a larger stellar radius can trade-off versus a larger impact 
parameter (i.e., transit at higher latitude) to produce the 
same transit duration. In the limit of a small planet and 
zero limb darkening, a chord at high latitude across a large 
star would produce the same transit curve as a chord at 
low latitude across a small star. This mild degeneracy 
is broken by the duration of ingress and egress, so that 
stellar radii exceeding approximately 0.23 Rq lie outside 
contours of acceptable x 2 . We derive a value for the plan- 
etary radius (R p ) by examining the R p values at points 
lying within the 1 a contour (that first contour is not illus- 
trated on Figure 3). This yields R v — 2.61 where 
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the asymmetric error range reflects the fact that the radii 
are not distributed symmetrically within the contour. We 
also derived best- fit values for the stellar radius and impact 
parameter, and these are listed in Table 1 for comparison 
to the Charbonneau et al. (2009) results. Our error range 
for the impact parameter (a factor of two) is relatively 
large, reflecting the degeneracy discussed above. However, 
our radius for the planet has a much smaller fractional 
variation within the Iff contour, and is in excellent agree- 
ment with Charbonneau et al. (2009). We conclude: 1) 
our value for the planetary radius is relatively insensitive 
to the mild degeneracy between stellar radius and impact 
parameter, and 2) the discovery value for the radius of 
this super-Earth (Charbonneau et al. 2009) is unlikely to 
be significantly affected by uncertainty in the optical limb 
darkening, or affected by stellar activity. 
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Fig. 1. Photometry of GJ 1211 during 5 transits using three different telescope/ wavelength combinations. Nominal transit curves are 
'•erlaid, based on the Charbonneau et al. (2009) system parameters, but shifted to fit the best central transit time. 
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FIG. 2. Top Panel: Raw J-baad photometry of the GJ 1214 transit observed from Kitt Peak using FLAMINGOS on May 29 UT. 1 he 
red line is a baseline fit using a 4th order polynomial, necessary due to the different effective wavelengths of the target and comparison 
stars when averaged over the filter bandpass. Bottom. Panel: Photometry divided by the baseline, with a theoretical transit curve (blue line) 
overlaid. The transit curve is based on IR limb darkening for M-dwarf stars {T — 3000/C fe logp -- 5.0) from Claret (1998), and using the 
Charbonneau et al. (2009) values for the system parameters. It was shifted in time to give the best fit, but, was not altered in depth or shape. 
These J-band data are publically available in the NStED data base. 
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fig. 3. Contours of constant x 2 for the best fit, to our J-band transit, in the plane defined by the stellar radius and impact parameter. 
The contours have been converted from x 2 to statistical significance, and are plotted for 2, 3, 4, and 5a significance. The cross shows the 
values and ±l«r range from Charbonneau et al. (2000), and the square indicates the position of our x 1 minimum. Note the substantial range 
within the contours where the same goodness of fit is obtained by trading a smaller stellar' radius for a lower crossing latitude (i.e., smaller 
impact parameter). This is due to a mild degeneracy in the transit fit, as discussed in Sec. 4. 
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Table 1 

Best-fit values for planetary and stellar radii, impact parameter, and bphbmerjs parameters, compared to the 

RESULTS OF CHARBONNEAU ET AL. (2009). 

~ParameteF BesTfittedl^ Clmrlx^^ 

~Planetary radius 2.6li^fji?a 

Stellar radius 0.20ll|;|gjR & 0.2110 ± 0.0097 R s , 

Impact parameter 0.175 + J'{*'J 0.354.to'og2 

Central transit time 2454964.94390 ± 0.00006 (UTC-BJD) 2454983.9087558 dfc 0.0000901 

Orbital period 1.5804043 ± 0-0000005 days 1.5803925 ± 0.0000117 

Note: Charbonneau et al. {2009} list only individual measured transit times, not an averaged epoch. We quote their most 
precise individual transit. 


Table 2 

Times of GJ 1214b transit center, given as UTC-based BJD. 


UT Date 

Telescope 

Filter 

BJD k Error (la) 

April 29. 2010 

UDEM 

I c 

2455315.79343 ± 0.00042 

Mav 29, 2010 

KP 2- m 

.1 

2455345.82126 ± 0.00014 

Mav 29, 2010 

KP VCT 

B 

2455345.82133 ± 0.00037 

June 6, 2010 

UDEM 

l 

2455353.72332 ± 0.00036 

June 17, 210 

UDEM 

I c 

2455364.78669 ± 0.00029 


Note: Adding these transit times to those of Charbonneau et al (2009), gives the best-fit ephemeris parameters (period 
and transit time) listed in Table 1. 


